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Enhanced Casimir effect for doped graphene 
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We analyze the Casimir interaction of doped graphene. To this end we derive a simple expression 
for the finite temperature polarization tensor with a chemical potential. It is found that doping 
leads to a strong enhancement of the Casimir force reaching almost 60% in quite realistic situations. 
This result should be important for planning and interpreting the Casimir measurements, especially 
taking into account that the Casimir interaction of undoped graphene is rather weak. 


Introduction Graphene, which is a two-dimensional 
sheet of carbon atoms possesses many unusual proper¬ 
ties and attracts a lot of attention. Particular excite¬ 
ment among the theoreticians is caused by the fact that 
the spectrum of quasi-particles in graphene is described 
by the quasi-relativistic Dirac model with the effective 
propagation speed of about 300 times less than the speed 
of light. This continuous model turned out to be very 
successful in describing a broad range of effects [1], for 
instance optical properties of graphene as the absorption 
of light [2] and the (giant) Faraday effect [3], to mention 
a few. 

In the recent years, the Casimir effect for pristine 
graphene was studied both for zero [4, 5] and finite [6- 
8] temperatures. For not too large temperatures (as 
compared with inverse distance between the interaction 
sheets) the effect between graphene monolayer and ideal 
metal is defined by the fine structure constant a ~ 1/137 
and is roughly 2.5% of the one between two ideal metal 
plates. Such small forces are on the limit of sensitivity of 
modern experimental techniques. For high temperatures 
(or separations) the effect is hugely reinforced [8], but the 
measurements under such conditions is a separate quite 
challenging task, which is not completely solved yet even 
for metals, see e.g. [9]. It is not surprising therefore that 
just a single experiment has been performed until now 
[10]. This experiment revealed [11] a good agreement 
with the theory [12]. Possibilities, opened by doping, 
were however not explored there. 

Previously the Casimir interaction of doped graphene 
was studied in [13, 14]. The reflection coefficients used in 
that works were expressed though quantities whose ex¬ 
plicit dependence on the temperature remained unknown 
and the results of Refs. [13, 14] are mutually contradict¬ 
ing. Therefore, specifically after experimental confirma¬ 
tion of the Dirac approach to Casimir energy of undoped 
graphene [10, 11], it seems to be important to extend the 
approach to doped graphene. 

In this letter we consider the Casimir effect at fi¬ 
nite temperature, chemical potential and mass, and find 
a substantial enhancement of the effect in graphene- 
metal systems which potentially permits to avoid the 
above mentioned experimental difficulties. Our findings 
show that for relatively highly (but still feasibly) doped 


graphene monolayers the effect gets stronger by approxi¬ 
mately 60%. In the (formal) limit of an infinite chemical 
potential the Casimir interaction becomes 1/2 of that for 
the ideal metal. Our calculation is based on a complete 
representation of the polarization tensor of the fermionic 
quasi-particles in graphene at finite temperature, chem¬ 
ical potential and mass gap. The obtained result is sur¬ 
prisingly simple and can be easily analytically contin¬ 
ued to whole complex frequency plane, including, impor¬ 
tantly, real optical frequencies. It is based on the QED 
formalism applied to graphene-like systems in [15, 16], 
and generalizes the results of [4, 8, 17] and those of [18], 
which also allows continuation to the whole plane of com¬ 
plex frequencies, to the case of simultaneous presence of 
finite chemical potential, finite temperature and non-zero 
mass-gap. We also describe a very precise approximation 
scheme that considerably simplifies the computation of 
Casimir interaction with graphene. 

Model The theoretical description of the electronic 
properties of graphene based on the continuous Dirac 
model with a 2 + 1-dimensional action. In the notations 
of [19] it reads 

Sb = J d 3 x‘ij}(T l (idi - eA t ) - m)ip, ( 1 ) 

where l = 0,1,2 and x = (a: 0 ,a: 1 , a; 2 ). The gamma ma¬ 
trices J 1 are rescaled, 7 0 = 7 0 , 7 1 ’ 2 = vp”/ 1 ’ 2 , 7 q = 
— ( 7 1 ) 2 = — (7 s ) 2 = 1. We use natural units h = c = 
kfj = 1, and the Fermi velocity is vf ~ (300) -1 . Thereby 
we assume that the graphene monolayer is placed at the 
(a; 1 , a; 2 ) plane. The electromagnetic potential A M is nor¬ 
malized in such a way that e 2 = 47 tq: = ^. 

Reflecting the spin and valley degeneracy in graphene, 
the gamma matrices, 7*, are 8x8 being a direct sum 
of four 2x2 representations (with two copies of each of 
the two inequivalent ones). The value of the mass gap 
parameter in and mechanisms of its generation are under 
discussion [15, 16]. 

As shown in many previous works, see reference in [19], 
the electronic properties of graphene in the formalism of 
Quantum Field Theory (QFT) can be described by the 
one loop polarization operator, which is gauge invariant 
in QED-like model defined by (1). In Minkowski momen- 
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turn space it is given by 

H m " (p) =ie 2 J ^ tr (s(q)^ m S(q - p)T ) , ( 2 ) 


where p = (j>o,Pi,P 2 ), 9 = (qo, 9i, 92 ), and <5 is the causal 
(Feynman) propagator of the quasiparticles in graphene, 


S(qo,q.\\) = - 


(<7o + h)lo - v F ^ - m 
(Qo + P + mo ) 2 - v 2 f q/ - m 2 


(e > 0). Note that due to the quasi-relativistic nature 
of excitations in graphene S also depends on the Fermi 
velocity v F . Further notations are q B = (g 1 ,^ 2 ), cfy = 
7 1( ?i + 7 2 < 72 , and p is the chemical potential. 

Temperature is introduced using the Matsubara for¬ 
malism. In the 7 -trace in (2), which can be calculated 
immediately, see e.g. Eq. (A20) in [17], one has to sub¬ 
stitute the integral by a sum, 


-2nT J2 > 


27riT(fc + 1/2), 


where k is integer. The external frequency of the po¬ 
larization operator is bosonic, po —> ip 4 = 2iriTn, 
n = 0,1,2,... 

Calculation of the polarization operator for finite tem¬ 
perature, mass gap and chemical potential All compo¬ 
nents of the polarization tensor can be expressed via two 
scalar quantities (form factors), n tr = n 0 o — Iln — II 22 
and n 00 [ 8 , 20]. As in the QED/QCD cases these quan¬ 
tities consist of the vacuum part and a part carrying the 
dependence on T and p. 

n xx ( P;/ x,T) = n^ ac) (p) + AU xx (p;p,T), (5) 


Here the distribution function, 2 = (el >ir+ 'd/ T + l ) -1 + 
( e (j<-/i)/T _|_ i)- 1 ^ carries the dependence on T and p. 
Further notations in ( 6 ) are 

M 0 o = -P 2 + ^iPiX + 4>r 2 , 

M tr = -p 2 + 4>r(l - Vp)(ip 4 + x) + 4 v F m 2 , 

Q=p 2 ~ 2i P 4X, p 2 =pl+v 2 F p 2 { , Ph = |Ph|. 

Note that An xx , ( 6 ), does not have UV singularities. 

For the vacuum part, n( vac ), one can directly use the 
well-known expressions [4] valid for graphene, 

j-r(vac) __ a ^P\\ Tr(vac) _ at<P{p~ + p 2 ) ^ 

11 00 - p 2 ’ 1A tr - p 2 > \<) 

where $ = 4 m + ^ ~ ~ m arctan • From now on 

we set m = 0 (gapless graphene) unless otherwise stated. 

One of the advantages of the decomposition (5) is the 
absence of the summation over Matsubara frequencies, 
which permits relatively easy derivation of the limiting 
cases. In particular, in the limit of zero temperature (but 
not zero chemical potential) we obtain 


An, r = 8 a x 


p 2 +p 2 


\v F L 4 p 
(1 - v 2 F )p\ 
4 p 

AHnn = 8 a x 


^ — ^4lm xy/x 2 + 1 + l°g (x + \/x 2 + 1 


[x + v x 2 + 1 


xy x 2 + 1 


v F 4 p 


( 8 ) 


(9) 


where xx stands either for ’tr’ or ’00’. The vacuum part, 
n xx lC ' 1 (p) , corresponds to p = T = 0. While such decom¬ 
position is a well known feature of polarization tensor in 
different theories, see e.g. [ 21 ], its realization in partic¬ 
ular cases and derivation of simple transparent formu¬ 
las may be a challenging task. One transforms the sum 
over the Matsubara frequencies (4) into a contour integral 
consisting of three parts, one of which gives the original 
integral over the continuous (Euclidean) momenta and 
the other two can be taken explicitly by the Cauchy theo¬ 
rem. The remaining integral over the in-plane momenta, 
Qii = (?i, 92)1 can be further simplified by performing the 
angular integration. Recently such procedure was ap¬ 
plied to graphene at p = 0, T ^ 0 in [18] and T = 0, 
p 0 in [17]. Omitting the technicalities we arrive at 

An xx = (6) 



with x = 2lfl Pi . The formulae above are the analogue 
of the (B2) [17] taken at Matsubara frequencies, and one 
can check that in the appropriate limits they reproduce 
the results of other authors [4, 8 , 15-18]. Similar to the 
results of [18], the representation (5) with ( 6 ) directly 
permits for continuation to the real frequencies, and thus 
can be applied for investigation of the optical properties, 
surface plasmons and other effects in graphene for finite 
temperature and chemical potential. 

Enhancement of the Casimir effect The Casimir en¬ 
ergy density (per unit area) for two parallel interfaces 
separated by the distance a is given by the Lifshitz for¬ 
mula [ 22 ] in terms of the reflection coefficients, T’xe TM , 

r TE tm> °f the TE and TM electromagnetic modes on the 
two interfaces, 

OO « 72 

« = k ° T E / WT E Ml - e- 2 '-r«r< 2 >), 

n——oo J X—TE,TM 

( 10 ) 
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where p = ^p\ + p 2 , and P 4 = 2ixnT are the Matsubara 

frequencies 1 . The reflection coefficients are taken at Eu¬ 
clidean momenta r = J"(p 4 ,Pu). They were derived in [8] 
in terms of the polarization operator components n 0 o,trj 

_ piioo _ p 2 n 00 - p] n tr 

r ™ piioo + 2 rTE p 2 n 00 -pif(n tr + 2p)’ 

(11) 

(and rederived in numerous papers afterwards). For the 
perfect conductor, rv^ = 1, = —1. Combining (10) 

with (11) and using (5,6) for the polarization operator 
at finite temperature and chemical potential, we are able 
to calculate the Casimir energy density, £, the Casimir 
pressure T = — d£/da , and its’ gradient, Q = dT/da = 
—d 2 £/da 2 , between a doped graphene layer and an ideal 
metal plate. 

It is instructive to consider the case of very large p 
first. In the formal limit p —» oo, both n 0 o and II tr have 
identical asymptotics, 

8 a 

n xx cs ° 4 p + ..., (12) 

/J->oo Vp 

which can be readily deduced from (8,9). Thus the elec¬ 
tronic properties of graphene would be expected to be¬ 
come closer to those of ideal metal. However, due to the 
specific structure of the reflection coefficients (11), it is 
only the contribution of TM mode to the Casimir inter¬ 
action which grows in this limit. Thus, at p —> oo, the 
Casimir interaction reaches the value of 1/2 of the ideal 
metal - ideal metal one, 


£ = 
fl—tOO 



k B T C(3) 
167TC1 2 


(13) 


The same result was obtained for the high tempera¬ 
ture limit [8]. Interestingly, the Casimir interaction of 
graphene can never reach 100% of that for the ideal 
metal. Eq. (13) gives a very rough idea on how far the 
enhancement of the Casimir effect with p might go. Prac¬ 
tically, it hardly makes sense to consider \p\ exceeding a 
couple of eV in the framework of the Dirac model. 

Our numerical analysis show that at distances about 
100-300 nanometers the Casimir effect between a per¬ 
fect metal plate and doped graphene is highly enhanced 
even for relatively moderate values of the chemical po¬ 
tential. On Fig. 1 we compare the behavior of the ratios 
of the energy, pressure and the pressure gradient at given 
values of the chemical potential to the corresponding val¬ 
ues for pristine graphene, as functions of distance. For 
p = O.leV, the interaction force for the case of doped 
graphene is only 3.5% higher then in the case of pristine 


1 To restore physical units in Eq. (10) it is enough to apply the 

conversion rule (eV)" 1 ~ 1.97 ■ 10 _7 m. 



FIG. 1. The ratios of Casimir energy density (full lines), 
pressure (dashed lines) and its’ gradient (dotted lines) at 
p = 0.1, 0.5, 0.8eV (black, blue and red lines, correspondingly, 
colour online) and at p = 0, between a perfect metal plate and 
doped graphene, as a function of distance a, [nm]. 


graphene. However, already for p = 0.5eV their ratio 
has a maximum of 32% at approximately 70nm, and for 
p = 0.8eV the force between a doped graphene layer and 
an ideal metal is almost 60% higher then that for a pris¬ 
tine one. All our simulations are performed at T = 300K 
and m = 0. 

One further notices, that the effect is more pronounced 
the more derivatives we calculate of the energy. Thus, the 
ratio of the energy at p = 0.8eV to its pristine value is 
1.52 at maximum, for the force T it is 1.54 and for the 
gradient Q it reaches 1.56. Moreover, at larger distances 
(400-1000 nm), the enhancement effect for the gradient 
diminishes much slower then that for the energy, which 
suggests preference for gradient force experiments. 

Reaching the values of the chemical potential of 0.8eV 
and higher might be a challenging task and would re¬ 
quire preparation of special samples. Without special 
treatment, the chemical potential in epitaxial graphene 
layers stays low up to the level of 0.3-0.4 eV or smaller 
[23, 24], However, the Fermi energy shifts of order 0.8 
eV are achievable in epitaxial single layer graphene due 
to molecular doping [25]. Under certain circumstances, 
doping may lead to considerable inhomogeneities in the 
charge distribution, see e.g. [26], that may give rise to 
additional forces in the Casimir experiments. Due to the 
strong charge density dependence on the nature of ac¬ 
ceptor/donor mechanism, these forces should be treated 
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FIG. 2. The ratios of Casimir energies between a perfect 
metal plate and graphene, as a function of distance a, [nm], 
m = O.leV, fj, = 0.8eV. 


individually for each particular experiment. 

The distance dependence of the Casimir energy for 
doped graphene will be altered as compared to the pris¬ 
tine one as it is usual when an additional dimensionfull 
parameter is introduced. However, its detailed study is 
beyond the scope of the present letter. 

Influence of the mass gap. From the physical point 
of view, the larger is the mass parameter in (1) the less 
should be the conductivity of quasi-particles and, conse¬ 
quently, the smaller the Casimir effect. This was shown 
explicitly in [4] for T = y = 0. One can also show that 
the influence of mass is negligible as far as m <C y. In par¬ 
ticular, in the formal limit (12) any dependence on the 
mass disappears. Our numerical simulation shows, see 
Fig 2, that for y = 0.8eV and m = O.leV doping gives 
up to 70% enhancement of the Casimir energy density 
(red line, to be compared with the full red line of Fig. 1). 
In the same picture it is shown that the influence of the 
mass gap on the value of the energy for doped graphene 
(blue dotted line) is almost negligible, while the energy 
for pristine graphene gets lower by about 15% (blue full 
line). Therefore, doping becomes even more important 
for gaped graphene. 

Approximating the Casimir energy. For finite tem¬ 
perature, the numerical calculation of the Casimir energy 
in a realistic set-up requires summing a large number of 
contribution to the sum over the Matsubara frequencies. 
Following the ideas of [8] (which were later confirmed in 
[27]), one might facilitate significantly this calculus by 
considering the T = 0 approximation (8,9) for the po¬ 
larization operator in all terms of the summation (10) 
except in the zeroth one. The comparison of exact re¬ 
sults for the Casimir energy £ with such approximation, 
£ app , is given on Fig. 3. As one can see, the error is 
less then 0, 5% for the pristine graphene and one order 
of magnitude smaller for doped one. This confirms once 
again at an even better level the asymptotic considera¬ 
tions delivered in [8]. 



0 50 100 150 200 250 a [ nm ] 300 

FIG. 3. The ratios of Casimir energy to its approximation 
according to [8], £/£ apP , for y = 0, 0.5, 0.8eV (full, dashed 
and dotted lines respectively) between a perfect metal plate 
and doped graphene, as a function of distance a, [nm]. 


Summary. In this letter we calculated the polariza¬ 
tion operator for the quasi-particles in graphene at non¬ 
zero temperature, chemical potential and mass gap ap¬ 
plicable at all complex frequencies without a need for any 
special procedure of analytical continuation. This result 
can be used in a variety of physical problems, including 
investigation of TE surface plasmons in graphene [17], 
quantum reflection [28], Casimir interaction, etc. 

Basing on these results and the Lifshitz formula, we 
numerically simulate the Casimir interaction between a 
doped graphene monolayer and an ideal metal. For high 
but feasible doping we predict the enhancement of the 
Casimir effect as compared to the case of the pristine 
graphene of up to 54% for the Casimir force, and of 
up to 56% — for the force gradient. The high doping 
of graphene is shown to bring significant enhancement 
in the values of the force gradient at a wide range of 
separations which should facilitate future experimental 
measurements. At such levels of doping the influence of 
the mass-gap is not important. We saw also that even 
moderate values of the chemical potential have a non- 
negligible effect on the Casimir force and thus should be 
taken into account in realistic description of experiments 
together with real material properties, finite temperature 
and mass-gap parameter, if present. 

All calculations of the present paper were performed in 
fully retarded approach valid at all distance. It may be 
interesting to study whether the non-retarded approach 
may deliver a good approximation at some distances. 

Finally we note, that the considerations given in this 
letter are concerned with the proper graphene proper¬ 
ties and the enhancement of the Casimir interaction is 
invoked by the change in its conductivity. Thus, we can 
conclude that even in the experiments involving a real 
metal and/or graphene on a substrate the enhancement 
effect must be present, though its particular value may 
differ from the one presented here. The good concor¬ 
dance between the force gradient measurements [10] and 
theoretical considerations presented in [11] shows that 
the graphene samples used in [10] were rather pristine. 
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